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1. Introduction 



Two loop calculations are not a rarity anymore in chiral perturbation theory ^, 
|, |, 1^, |, |, 0, |n|, 0, 0, |14|, |16|, 0, 0, |0|, |23, |2|. Due to the presence 
of masses the calculations are rather complex and only few explicit analytic results are 
available, particularly in the three flavour case. 

In cases where the final result depends on well separated mass scales it is reasonable 
to expect that a few terms in the expansion in terms of the ratios of the scales yield a good 
approximation to the full result. In turn, there is hope that the individual terms in the 
expansion represent simpler expressions and may even be given in closed form. Because 
of the complexity of the full results, it is desirable to have a method which allows an 
expansion prior to the integration. Moreover, integral representations suitable for numerical 
integration for generic momenta sometimes exhibit singularities at the (pseudo) thresholds. 

Below, we discuss integrals that contribute to the ttK scattering amplitude at threshold 
and thus depend only on the pion and kaon masses M,r and Mk (throughout, we work 
in the isospin limit m„ = md)- The envisaged physical application is to find analytic 
representations for the expansion of the ttK scattering lengths in M-,^/Mk- The most 
precise values for the ttK scattering lengths are presently obtained from an analysis of 
Roy-Steiner equations |24, 25|. Alternatively, particular combinations of ttK scattering 



lengths may be extracted from experiments on ttK atoms [26, 27, 28|. The nK atom 
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decays due to the strong interactions into ir^K^ and a lifetime measurement will allow one 

1 1/2 3/2 

to determine the isospin odd S-wave ttK scattering length Oq = ^(ag — Oq )• Such a 
measurement is planned at CERN [^^. The scattering length Cq differs from other low- 
energy parameters in irK scattering in that it is protected against contributions of in 
the chiral expansion: In the framework of SU(2) chiral perturbation theory |3C, |32| , 



where the kaon is treated as heavy, there exists a low-energy theorem |^2| which states 
that the Weinberg current algebra result [33, 34] receives corrections of order only, 



{1 + 0(M^^)}. (1.1) 



" 87rF2(M^ + Mi^) 

Here M^-, Mk and i^V denote the physical meson masses and the physical pion decay 



constant. It was therefore expected that the one-loop result |35, 3C, 37] in the SU(3) theory 
]|38|] should represent a decent estimate for this scattering length. However, the dispersive 
analysis from Roy-Steiner equations ]25] and the chiral two- loop calculation ]|l^ are not 



in agreement with this expectation. In fact, the numerical analysis performed in Ref. ]18] 
showed that the two-loop order corrections to Oq are of the same order of magnitude as the 
one- loop contributions. To understand the nature of these rather substantial corrections, 
a partial analysis of Oq at next-to-next-to-leading order in SU(3) chiral perturbation was 



performed in Ref. |3£]. The present work represents a first step towards an analytic 
two- loop representation for the expansion of the irK scattering length in M-j^/Mk- 

The technology we employ to obtain our results is known by the name of expansion by 
regions ]^ , ^\\ . The recipe for the method is simple: i) In the integrand identify all the 
relevant integration regions, ii) expand the integrand in each region, iii) sum up the integrals 
(integrated over the full integration domain) of all the contributions. The method obviously 
owes its applicability to the fundamental property of dimensional regularization, being that 
the integral only receives contributions from scales which are present in the integrand. The 
divergences encountered in intermediate steps are also tamed by dimensional regularization. 
The number and nature of the regions to be considered depends on the problem at hand. 
In the present case, it is sufficient to consider hard (Mk) and soft (Mjr) regions. Since 
we furthermore do not violate manifest Lorentz invariance we only need to consider two 
distinct regions for the loop momenta. 

In Section ^, we introduce the method for the example of the two point integral J 
which contributes to irK scattering at the one loop level. On the basis of the comparison 
with the expansion of the known full result we are able to demonstrate the validity of the 
method in this case. The prescription translates to the case of two loop integrals which are 
considered in Section ^. Specifically, we apply the method to three variants of fish type 
(two loop vertex) integrals. We demonstrate the usefulness of the method by computing 
the two leading terms in the expansion of m/M of each integral. Apart from products of 
one loop integrals (see Appendix |A|), the determination of the coefficients in the expansion 
involves several variants of on-shell sunset type integrals whose evaluation is described in 
Appendix]^. Finally, Section ^ contains our conclusions. 

In the present work, we are making use of methods that are common in fields other 
than chiral perturbation theory. While we are unable to give a comprehensive bibliography 
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on the subject here, this is mitigated by the existence of excehent reviews |41, 42 1. Ref. |40| 
is equahy suited for an introduction to the subject of the expansion by regions. Note that 
while at present there are no proofs for the vahdity of the method in the Minkowskian 
region, there are also no known counterexamples. The expansion in a small mass ratio 
m/M was studied in Ref. in the case of the two loop two point function and in the 
framework of nonrelativistic QED in Refs. |44, 45], with techniques very similar to ours. 
For two loop integrals with arbitrary masses, see Refs. 48, 49 1. Last but not 

least let us mention the existence of other methods suitable for the expansion of Feynman 
integrals in small parameters. A prominent example is the Mellin-Barnes representation 



llg, p, n, m, g n, n, 113, |2|. 



2. J Integral 

As an illustration of the method, we consider the J Integral at threshold^ {p + P)^ 

(m + M)2, 

d'^k 1 1 



Jthr(m,M) 



{2-kY 

Jthr(0,M) 



-2pk - 
1 



A;2 - iO 2Pk - A;2 

,2 



(2.1) 



m 



■In. 



m 



+ 0((i-4) 



(47r)2 m + M 

In order to simplify our notation, we use throughout the symbols M = Mk and m = M^. 
As shown in Fig. |l|, we let the external kaon momentum = flow through the kaon 
line so that the kaon mass is cancelled. Otherwise the momentum assignment is free. We 
choose a symmetric configuration where p'^ = re? denotes the external pion momentum. 
Further our reference frame is chosen such that P ^ M and p m (center of mass system) . 
The regions for the loop momenta k are defined as follows 



hard (h) : 
soft (s) : 



k M , 
k ^ m . 



(2.2) 



The expansion of the hard region {k ~ M) looks as follows: 



i 



d'^k 



{27tY -k^ 

oo 



iO 2Pk -k'^-iO 



E 

n=0 



2pk 



-/c2 - iO 



(2.3) 



n=0 



Taking the measure into account, the resulting order, = 0(m"/M"), is anticipated. 
Studying the behaviour at zero and infinity shows that the leading order term Jq , coinciding 
with Jthr(0, M), contains an ultraviolet divergence and is infrared finite. Explicitly, we get 



■^0 



Jthr(0,M) 



-2e 



(47r)^ 



M2 

ln-^ + 2)+0(e) . 



(2.4) 



^An analogous consideration for the on-shell integral, [p + P)^ = M^, may be found in Ref. 
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P-k 




p + k 



Figure 1: J Integral: The solid line denotes the kaon, the dashed line stands for the pion. 



We indicate the nature of the divergence by e — > euV) where 

e = 2--, /z2 = ^_/z^ (2.5) 

For n > 1, the hard part is ultraviolet finite and infrared divergent. The next-to-leading 
order takes the form 

^N-^i^f-'l^^ + ^UoW, (2,6) 



Af(4T)2VeiR 

where e — > eiR. 

The expansion of the soft region (A; ~ m) takes the form 

_ 1 r 1 l__ ^ 

oo 



n-l 



i2-KY -2pk -k^ -iO 2Pk - iO ^ \ 2Pk - iO 

= T.Jn- (2.7) 

n=l 

Again, = 0{rn"' /M^). In this case there is an ultraviolet divergence for all n(> 0) and 
the integrals are infrared finite. For the leading term in m/M we find 

This result is obtained easily by choosing the frame such that P = M/mp and subsequent 
partial fractioning. In the sum of the hard and soft part at order m/M 

the infrared pole term in and the ultraviolet pole term in have cancelled (since 
euv = eiR = e) to give the desired result as in Eq. (^). The cancellation of an infrared 
with an ultraviolet divergence is in fact common in dimensional regularization, as the 



-4- 



following example shows |55]: 
1 



1 



+ 



m 



fc U-A;2 -i0)(m2 - A;2 -iO) (-A;2 - i0)2(m2 - P - iO) 
1/1 1 
euv eiR 



(4vr) 



. 



(2.10) 



In the following, we will refrain from detailing the nature of the divergences encountered 
at intermediate steps and trust in dimensional regularization to take care of the necessary 
cancellations. This procedure works on the premise that the divergence structure of the 
end result is known from other sources. This is of course the case for the J integral as well 
as for the two loop diagrams considered in the next section. 

To conclude the present section we note that it is straightforward to establish the 
relations 



J, 



n+l 



M " ' 



m 



M " 



(2.11) 



hereby extending the validity of the approach in the case of the integral Jthr("i, M) to all 
orders in m/M. The method is equally valid when the dimension d is left arbitrary. 

3. Fish type integrals 

In this section, we discuss the two loop fish type integrals occurring in ttK scattering at 
threshold, (p + i-*)2 = (m + iVf)2. Our aim is an expansion of the fish type diagrams in 
powers of the mass ratio m/M. To achieve this we perform an expansion by regions. The 
set of regions that contribute to the expansion of the fish type diagrams at threshold in 
general consists of 5 regions [45|, 



h-h : 


k 




M , 




k- 


l ~ M 


h-h' : 


k 




M, 


I M , 


k- 


I ~ m 


s-h : 


k 




m , 


I M , 






h-s : 


k 




M, 


I ^ m , 






s-s : 


k 




m , 


I ^ m . 







(3.1) 

As indicated, the difference in the h-h and h-h' regions lies in the fact that the difference 
k — lis counted as hard and soft, respectively. Note that in the examples considered below, 
the h-h' region does not contribute. The loop momenta k and / are shown for example in 
Fig. |2|. In the following, we discuss this topology with different particles running in the 
subgraph. 

3.1 Pion loop 

Let us start with the fish diagram containing an internal pion loop (for our notation, see 
Appendix |^ 

1111 



Fl 



kim? -P -iQm? - {k- If - iO 2Pk - k'^ - iO -2pk - k'^ - iO ' 



(3.2) 
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p + p 



p ■■ 



p + k 



Figure 2: Fish diagram with an internal pion loop. 



shown in Fig. g. The fish type diagram contains a nonlocal divergence, which is known 
for arbitrary masses and external momenta, see Eq. ( |B.16| ). At threshold no additional 
infrared singularities occur. 

The h-h region contributes at order one, the s-s part is of order m/M, while the 
h-h' and h-s regions start contributing at order rm? /M"^ which is beyond the scope of our 
calculation. The expansion in the s-h region leads to scaleless integrals so that this part 
vanishes altogether. We thus have 

Fl = Fl^i^ + Flf + Flf + O f ^) . (3.3) 



The expansion of the h-h region {k ~ M, I ~ M) looks as follows 

11 11 

-P - iO -{k - ly - iO 2Pk -k'^ -iO -k"^ - iO 



Fl 



hh 



1 + 



2pk 



-fe2 - iO 



+ 



m 



Fl"^^ + Flf + O 



m 



(3.4) 



To evaluate the h-h contribution at leading order, we first carry out the integration over 
the subgraph momentum I. The resulting integral over k is then again of the form of a one 
loop integral and the result can be written in terms of one loop functions. 



Fl^^ 



M' 



2d-8 



/i,i{0,0;-l}/3_d J0,1;1} , 



(3.5) 



where the Ia,f3 are defined in Appendix^. At leading order in m/M the fish type integral 
Fl is given by 



Fin = Fl 



hh 




(3.6) 



-4e 



1 



1 

(47r)4 \2e^^Ye 



2 In. 



AP 

^2 



M2 
/x2 



M2 
/x2 



+ ln2 — - 51n — + y + — J> + 0(e) , 



57r2 
12 



where the parameters e and fi are defined in Eq. ( |2.5D . For the next-to- leading order 
contribution, we find 



Flf 



m 

'm 



Fl 



hh 



(3.7) 
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p + p 



p ■■ 



p + k 



Figure 3: Fish diagram with an internal kaon loop. 



To evaluate the contribution from the s-s region (A; ~ m, / ~ m) at order m/M 

111 



Fl 



ss^ f L 



(3.8) 



P - iO m2 - (/c - /)2 - iO 2Pk - iO -2pk - k^ - iO ' 
we exploit the fact that our reference frame can be chosen such that p = m/M P. Now we 



may use the following trick |61|] : In the integral we insert the identity 



1-iA, 

ddk,," 



(3.9) 



and integrate by parts to obtain a relation between the s-s part and sunset type diagrams 

,2d-8 



Flf 



m m 



M d-A 



[25i,i,2{l, 1, 1; 1} + 52,i|i{l, 1; 0, 1} - /2{l}/i|i{l; 1, 1}] • (3.10) 



The one loop integrals and Iq,|^ are listed in Appendix The sunset 52.i|i{l, 1; 0, 1}, 
given in Appendix]^, can be expressed in terms of Gamma functions. To evaluate the finite 
part of FVf we further need the equal mass case on-shell sunset S'i,i,2{l5 1, 1; 1} at order 
e. The evaluation of this sunset type diagram is outlined in Appendix ^ and the result at 
order e is listed in Eq. ( [8.81) . 

If we add the contributions from the h-h and s-s regions, the result for the fish diagram 
-Fl at order m/M yields 



Fh 



Flf + Flf 



m ^ 



-4e 



1 , M2 

- -21n — 
e 



M (47r)4 
3.2 Kaon loop 

The fish graph shown in Fig. ^ looks as follows 
F2 ' ' ' 



2 2 2 

, m , o m , m r, , ^, ^ 



(3.11) 



kiM^ -P -iOM^-ik- ly - iO 2Pk - k^ - iO -2pk - k^ - iO ' 



(3.12) 



The leading contribution stems again from the h-h region, the s-h part starts contributing 
at order m/M, while the h-h', h-s and s-s contributions all vanish identically. 
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The expansion of the h-h region {k ~ M, / ~ M) takes the form 

r 1 1 1 

IkiM"^ -P - iO M'^-{k- ly - iO 2Pk - k'^ - iO -k"^ - iO 



F2 



hh 



1 



X U + 



2pk 



F2^^ + F2T + O ( TTo 



-A:2 - iO 



hh 



+ 



/ 2 

/ m 



m 



VM2 

2 \ 



AP J 



(3.13) 



As with the s-s contributions to the pion loop diagram discussed above, we may use inte- 
gration by parts to reduce F2^^ to the sunset type diagram 5'i_i^2 



F2, 



hh 



1 



d-4 

At leading order F2 amounts to 



[25i,i,2{M, M, M;M^}- J2{M}/i,i{0, M; M^}] 



(3.14) 



F2o = F2^^ 



(3.15) 



-4e 



1 1 

+ 



(47r)4 [2e2 2e 



5 - 2 In ■ 



7? 



M2 



M2 19 Tvr^ 



+ ln^__51n_+ ^ 



12 



+ 0(e) . 



At order m/M, both the h-h and s-h regions (k m, I ^ M) contribute to F2. The h-h 
part takes the form 



F2f' 



m 
M 



i^2^'^-5i,i,2{M,M,0;0} 



(3.16) 



The result for the sunset integral 5'i^i.2{M, M, 0; 0} can be given in terms of Gamma 
functions and is listed in Appendix The s-h contribution yields 



F2f 



I 

Jkj 



m 



1 



In total, we get for F2 at next-to-leading order 



[M2 - /2 _ io]2 2Pk - iO -2pk -k'^-iO 

2 ^21 



/2{M}/i|i{m;m^m^} . 



F2i = F2f' + F2f 



m /i 



-4e 



1 M2 

- -21n^ 
e /x^ 



, m2 1 , m2 , m2 27r2 
'^AP^-2'^JP-''^JP-'^ — 



M (47r)4 
3.3 ttK loop 

The fish diagram with an internal ttK loop, shown in Fig. ^, is given by 

^ f 1 1 1 1 

I 2pl -P-iO -2P{k -l)-{k- If - iO 2pk - k^ - iO -2Pk - k'^ - iO 



(3.17) 
(3.18) 

+ 0(6). 



. (3.19) 



The leading order contribution is generated by the h-h region, the s-h part contributes at 
order m/M, while the contributions from the h-s and s-s regions are of higher order in 
m/M and the h-h' part vanishes identically. 
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p + p 



p ■■ 



Figure 4: Fish diagram with an internal nK loop. 



The expansion of the h-h region {k ~ Af, I ~ M) looks as follows 

1 1 1 



F3 



hh 



-2P( 



X 1 



2pk 



l)-{k- ly - iO -A:2 - iO -2Pk - k"^ - iO 
2pl 



-k^-iO -P - iO ^ [^M^ ) 



2 \ 



F3^^ + F3f + o(^ 



(3.20) 



To evaluate the h-h part at leading order, we again make use of integration by parts, in 

1 d 



this case by inserting 1 = 2W~^tJ-- Integrals of this type were also discussed in Refs. [45|. 
The result takes the form 



M 



2d- 8 



d-3 

and F3 reads at order one 



/2{l}/i,i{0, 1; 1} - /i,i{0, 0; -1}/2_| 2{0, 1; 1} 



(3.21) 



F3o = F3, 



hh 




-4e 



(47r)4 [2e2 2e 



5 - 2 In ■ 



/.2 



M2 
^2 



M2 19 7r2 
^ ' 2 4 



+ ln^^-51n^ + ^-^[>+0(e) . (3.22) 



At order m the function F3 receives contributions from the h-h 



2d-8 



- /2{l}/i,i{0, 1; 1} + /i,i{0, 0; -l}/2-|,2{0' 1; 1} 
+ 2/i,2{0, 0; -1}/3„|_,{0, 1; 1} + 52,i,i{l, 1, 0; 0} 

and s-h regions (/c ~ m, / ~ M) 



F3f 



1 



1 



1 



1 



m 
M 



-P - iO 2Pl -P-iO 2pk -k'^-iO -2Pk - iO 

2 ^21 



h^i{0,M;M'}h\^{m;m',m'} . 



(3.23) 



(3.24) 
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The result for F3i can be obtained by inserting ^f-A;^ 
result for F3i amounts to 



l-{lf,-k^) in F3^\ The net 



d 



F3i 



F3f + F2,f 



m ^ 



-it 



M (47r) 



1 , M2 
- -21n^ 
e /i^ 



(3.25) 



In — ^ + -ln2 — ^-4 In — ^-1 \ + 0(e) . 



4. Conclusions and an application 



As discussed before, the divergent part of the fish diagrams Fl, F2 and i^3 is known. 
When subtracting the poles in d — 4 in a scale independent manner (cf. Appendix the 
results obtained in the previous section may be written in the compact form, 



— 1 

Fl 



F2 



F3 



(47r)4 
1 

(4^ 
1 



19 57r2 



+ 



m 
12 



In' 



m 
M2 



5 In 



m 



M2 



+ 



111 
12 

^2 



+ 



m 
M 



2 M2 



2 In 



19 vr" m 



2 



m 
M2 



4 In 



m 



m 
M2 
2 



27r^ 



M2 



7 + 



^2' 

1 

3 



m2 \ 
^ / m2 



(4.1) 



The equations imply that our calculation has reproduced the correct divergent pieces. 
We are not aware of an independent determination of the finite parts. We note that for 
m/M = /Mj^+ ~ 0.28 a suppression of the subleading terms is clearly observed, except 
in the case of F2 where the leading term happens to be small. 

These results are obtained by expanding the integrands for soft (m) and hard (M) 
loop momenta and subsequently summing the integrals of all contributions, to the desired 
order in m/M. The coefficients of the expansion are then given by homogeneous Feynman 
integrals. We demonstrate in some detail how those can be evaluated using integration by 
parts identities. The most difficult integral encountered is an equal mass on-shell sunset 
type graph which is is evaluated by standard methods in Appendix |b|. In addition to the 
two loop results we demonstrate the validity of the expansion by regions in the case of a 
one loop example to all orders in m/M. 

For the method to work, it is necessary to route the hard external momentum through 
the heavy particle line. In the present case, the conservation of Strangeness guarantees 
that this can always be done. For the same reason the problem of ttK scattering in QCD 
can be attacked also in an effective theory where the kaons are treated as heavy from the 
start [^].^ In fact, the method of the expansion by regions and the formulation in terms 
of an effective theory are very closely related pO| , |4l| ] . 

The approach followed in Refs. |3^, |6^ and in the present work is a different one: 
Since the ttK scattering amplitude is already known to two loop order in SU(3) chiral 
perturbation theory fl^ one would like to make use of those results to obtain the scattering 
lengths. The integrals discussed in the present work are believed to be the most intricate 



^Similar are the cases of Heavy Quark Effective Tfieory 
Tfieory [||, ||, |66[ || . 



or (Heavy) Baryon Cliiral Perturbation 
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ones in what concerns the expansion in Af^/M^' so that we hope to have contributed a 
significant step towards this goal. 

The result for the ttK scattering amplitude |]l8| at threshold also involves diagrams 
with rj mesons running in the loops. Note that the mass of the 77 is not independent but 
can be expressed in terms of the pion and kaon masses by virtue of the Gell-Mann-Okubo 
relation, 



Mi 



Ml) + 0(m2 ^^k) ■ 



(4.2) 



Also in this case one therefore deals with a problem with two mass scales and the expansion 
by regions proceeds exactly as in the examples considered above. It should furthermore 
be obvious that the applicability of the method is not limited to scattering problems but 
that it could as well be applied to a similar analysis of masses and decay constants or form 
factors in SU(3) chiral perturbation theory. As a first application we have calculated the 
leading strange quark contributions to the pion mass 



Ml 



{ruu + md)Bo 

Cl 



C2 
„tree 



loop 



(4vrFo)^ 
-16(47r)2(L^-2L^ 



. tree , loop 
C2 -\- <-2 



rrisBQ 
(AttFo] 



+ 0{inu,md,ml 



2 ArrisBo 
— m 7^ — 

9 3/Li2 



64(47r)4 [{-C[, + Cl, + ^Cl,)Fi + ALl{Ll - 2L^)] , 



12 



+ 



73 
16 



38 
81 



32 ..380 



-In^ 
9 3 



TUsBq 16 r 



(4.3) 



21n|£5l 



2 
9 



4\/2 



Im 



Li2 



where Li2 is the dilogarithm, Li2(z) = X^^i ^ 01- We have made use of the repre- 
sentation for the pion mass of Ref. |jT^, for the definitions of the coupling constants see 
Refs. [38, ^] (Note that our Cl carry mass dimension —2). We have also introduced the 
following auxiliary symbols for combinations of the Lj, 



C\ = (4vr) 



'-2 

c 



(4vr) 



1 3 

29L^ - -LI + mil 



6L7 + ua 



+ — L3 + 2Ll + -Ll 
2 2^ 24 ^ 3^ 



4{Ll + L-J + LI 



3 - {AttY [8LI + 2(L^2 + L3 



ULl - 2LI + 12LI - 6L7 + 2LI 



(4.4) 



Taking the scale dependence of the chiral coupling constants into account one 
verifies that the coefficients ci and C2 are separately independent of the chiral renormal- 
ization scale ^ (Note that the product tUsBq is QCD scale independent). The relevance 



'For the S-channel, there are six additional diagrams containing r; mesons. 
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of the formula Eq. ( [4.3D lies in the fact that it fully specifies the dependence of the SU(2) 
coupling constant B on the strange quark mass at two loops, 




The term ci agrees with the result of Ref . [^] . For y/m^B^ = 484 MeV and ^ = 770 MeV 
the logarithm in ci amounts to about 0.14. Neglecting the terms proportional to the Lj 
and Ci in C2 one finds a value of about 7.6. The bulk of this contribution stems from 
the term —380/81 ln{msBQ/ fj?) ~ 4.4 while the contribution from the last line amounts to 
about 2.6 (Im[Li2 ^+3^^ ] ~ 1.001). With those values and Fq = 92 MeV as estimates, the 
expansion takes the rather peculiar form B/Bq = 1 + 0.02 + 0.23 + . . .. The presumably 
important impact of the contributions from the coupling constants Li and Ci remains to 
be investigated.^ We believe that further such studies could provide substantial insight on 
the convergence of the chiral perturbation series. 
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A. Table of one loop integrals 

The present Appendix lists results for a number of Feynman integrals that appear in the 
main text. Since the one loop integrals occur with various configurations of the indices, it 
is convenient to list them here in general form, leaving also the dimension d arbitrary. In 
cases of a singular configuration of indices, the relevant result is obtained via a limiting 
procedure. For a more complete list of this kind, see e.g. Refs. Throughout, we 

use the following shorthand notation 



*Note that the relevant combination Cig — C20 — 3C21 is Zweig rule violating. In the estimate of Ref. ||7i 
it therefore fails to receive a nonvanishing contribution. 




(A.l) 



We encounter various integrals of type one point or two point function, 




(A.2) 



(A.3) 
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In this notation the J integral at threshold is given by 

J,^,{m,M)=h,^{m,M;{m + Mf} = \ . (A.4) 

(47r)2 "-3 m + M 

For vanishing masses, Mi = M2 = 0, or in the case of the integral at threshold with one 
vanishing mass, Mi = 0, p'^ = M| = M^, one has 

laA^Ap'} = ^ . ' v(\v(R\vU R\ V-iO "°-^ , (A.5) 

(47r)2 T[a)T{li)T[d - a- f3) 



1 T{a + f3-f)T{d-2a-f3) 



A different class of integrals involves a denominator which is linear in the loop momentum: 
For p = one finds 

4|,{M;0,p'n = -^^^^^^^fT^^ , (A.8) 

(47r)2 ^-L 

while for = M^, = itp and integer (3 one has 
B. Two loop integrals 

In this appendix we collect various results for the sunset and fish integrals. We use the 
shorthand notation 

1 r d'^k r dH 
B.l Sunset type integrals 

The general integral of type sunset,^ shown in Fig. |5|, takes the form 

1 1 

(B.2) 



jk,i [Ml -P- iO]» [Ml -{k- If - i^f [Ml -{p- kf - iOp ' 

The following Feynman parameter representation for the integral is obtained by first 
performing the integration over the loop momentum / in Eq. ( |B.2| ) and subsequently com- 
bining the two remaining denominators |75|, 



^The integral also goes by the names of 'sunrise' or 'setting sun diagram' 
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Ml, I 



Ml, I 



p 


' M2, k-l ^ 


P 


— 


> 


—*■ — 



M3, p-k 



p 1 


' M2,k-l^ 






> 





M3, p - fc^ 



'Mi, p' -k 



P-P 

Figure 5: The sunset and fish topologies. 



The combination depends on Mi, M2 as well as on the Feynman parameter x 



Mt 



Ml ^ M| 



1 . (B.4) 

\ — X X 

In the following special cases the sunset integral may be expressed in terms of Gamma 
functions, 

1 r(a + /3 + 7-d) r(a + 7-f)r(/3 + 7-f)r(f -7) 



T{a)T{m-2) 



{AttY r(a + /3 + 27 - d) 
5,,/3,^{0, 0, M; M^} = /„,^{0, 0; -1} , 1 JO, 1; 1} [M^ - iO]'^""-^-^ 



(B.5) 

where the last result follows directly from Eq. ( |B.2| ). 

The evaluation of the diagrams in the main text does also require the result for the 
equal mass on-shell sunset integral with one denominator raised to the second power, viz. 
Si,i,2{M,M,M;M'^^. Furthermore, its contribution comes with a prefactor — 4) so 
that we need the expression including contributions of 0{d — 4). While the result in 



question is implicit in Ref. |74|, we pursued a more direct approach following Ref. |75|: In 
the Feynman parameter integral representation write (e = 2 — d/2) 

1 d . 



y 



-i+t 



e dy 



and integrate by parts to obtain 



5i,i,2{M, M, M- M^} = [M^ - iO] 



d 



-2e 



X <jl + 2e(l-y)— In 



r(26) 
e(4^)4-2e 

y 



dx [x{l — x)]~ 



+ (1-2/)^ 



y 



x{l — x) 



dyy" 



+ (1-2/)' 



(B.6) 



(B.7) 



-2e 



x{\ — x) 

Now the parametric integrals can be done by expanding the integrand around e = 0. We 
obtain 



5i,i,2{M,M,M;M2} 



2(47r) 



4-2e 
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We observe that this result is of equal simplicity as the on-shell sunset itself 77, 12, 74 1 



+ 0(e2). (B.c 



Finally, we also encounter a sunset diagram with one denominator which is linear in 
the loop momentum 

f 1 11 

5,,/3|^{Mi,M2; = J^^^ [^2 _ _ _ [Ml _ (/ _ ky - iO]/3 [-2p'fc - iOp ' 

(B.IO) 

For Ml = M2 = M and ]3 = we have 

mM-n«'n- 1 r(a + /? + 1 - d) r(a + 1 - f )r(a + 1 - f )r(i) 
V;3|7i^W,M,U,p i- (4^)dr(a + /? + 7-d) 2r(a)r(/3)r(7) 

X -iOj'^-^-^-afp'^ -iO]-i , (B.ll) 



which for a = (3 = 1 agrees with the result given in Ref. [43|. To obtain it, combine the 
first two denominators and perform the / integration. After this, the formula Eq. ( |A.8| ) can 
be applied. 

B.2 Fish type integrals 

The general integral of type fish, shown in Fig. |5|, is defined by 

i^a,/3,7,4Afl,M2,M3,M4;p^p'^(p-p')'} = (B.12) 

11 1 1 



k^i [Ml -p- io]° [M| -{k- If - iQf [M| - (p - ky - myi [Ml - ip' - kf - io]^ ■ 

In this notation the integrals Fl, F2 and F3 discussed in Sec. |3| are given by 

Fl = m, M, m; M'^,nf, (m + Mf} , 

F2 = Fi,i,i,i{M, M, M, m; M'^,m^, (m + Mf} , 
F3 = Fi,i,i,i{m, M, M, m; m^,M'^, (m + Mf} . (B.13) 

By combining the third and fourth denominators the diagram is lead back to the sunset 
0, 



i^a,/3,7,5{^l'^2,M3,M4;p2,p'2,(p-p')'} = (B.14) 

^!^^iyzz^-\l-zY-^S^,f,^,+s{Mi,M2,M,-,pl}, 

with 

M^ = zMI + {1 - z)Ml - z{l - z){p - p'f - iO , p, = zp+{l-z)p'. (B.15) 
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This representation is useful e.g., for the calculation of the (nonlocal) divergent part of the 
fish diagram with unit exponents, 



Fi,i,i,i{Mi,M2,M3,M4;p2,p'2,(p-p')'} = (B.16) 

r('2 - ^) A/f^'^-s rci - 

Mf-^^ hAMs, M4; {p - p'f} + 7 , 7, ,r + 0{{d - 4)°) , 
(47r)2 d-6 [Airf 



which, in the equal mass case, is seen to agree with the findings of Ref. |12|. The formula 
shows that the divergent part is in fact independent of the mass content of the subgraph. In 
particular, this implies that the divergent pieces of the integrals Fl, F2 and F3 discussed 
in Sec. |3| are equal. We can therefore define a quantity Fdiv, 



div 



(47r)4 Ud-4)2 ^_4 
2m 



5 - 2 In ^ 



,2 M2 
+ In 5 In 



1 , M2 
+ln^ 



A 



m + M 

d-Fdiv/d^ = 0(d — 4), such that the limiting expressions 



Fn= lim(Fn-Fdiv) , n=l,2,3, (B.18) 
are both finite and scale independent. 
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